Proximity effect and its enhancement by ferromagnetism in high-temperature 

superconductor-ferromagnet structures. 
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We consider a bi-layer consisting of a d— wave layerd superconductor and diffusive ferromagnet 
with a domain wall (DW). The c— axis in the superconductor and DW in the ferromagnet are 
assumed to be perpendicular to the interface. We demonstrate that in such a heterostructure the 
inhomogeneous exchange field enhances the proximity effect. It is shown that, whereas in the 
absence of the exchange field the d— wave condensate decays in the normal metal on the mean free 
path I, the superconductivity penetrates the ferromagnet along the DW over much larger distances. 
This happens because the presence of DW results in a generation of an odd frequency triplet s- 
wave component of the condensate. The phenomenon discovered here may help to explain a recent 
experiment on high temperature superconductor-ferromagnet bi-lyers. 

PACS numbers; 74.45.+C, 74.50.+r, 74.78.Fk 



Proximity effect (PE) in superconductor/non- 
superconductor heterostructures remains a popular topic 
of both experimental and theoretical study for several 
decades. It is well known that the superconducting 
condensate can penetrate a normal metal (N) over long 
distances ~ ^at = y/D/2TTT, where D is the diffusion 
coefficient in N and T is temperature 1]. In the limit of 
low temperatures T ^ the length ^jv can be very large 
and this is a long range PE. 

The situation is different in superconductor- 
ferromagnet (S/F) heterostructures. If the exchange 
field h in the ferromagnet is homogeneous, the Cooper 
pairs penetrate a diffusive ferromagnet over considerably 
shorter distances £_h = y/D/h. Usually T <C /i so that 
the length ^/j is much shorter than 0, 0] and PE is 
considerably suppressed with respect to the one in the 
superconductor-normal metal (S/N) structures. 

It turns out, however, that PE in S/F heterostructures 
can be long ranged provided the exchange field is inho- 
mogeneous. In this case a triplet odd frequency s-wave 
component of the superconducting condensate is gener- 
ated and it can penetrate into the ferromagnet over a dis- 
tance of the order ^tv, thus making PE similar to the one 
in the normal metals [3, 13| • As a typical inhomogeneity 
one can think of either thick 0] or thin ^] domain walls 
(DW) parallel (or perpendicular Q) to the interface or 
spin-active interface [5| . The prediction of the LRTC has 
been confirmed in several recent experiments ^] . 

The long range PE was usually studied in S/F 
heterostructures with conventional superconductors 
(s— wave singlet superconductivity). At the same time, 
S/F heterostructures can be fabricated using high tem- 
perature superconductors (HTS) [13, Hi, and study of 
PE in such systems is also important and interesting. 



In a recent work, the authors of Ref. [ll] using the 
scanning tunnelling spectroscopy measured the density of 
states (DOS) i'{e) at the outer surface of the ferromagnet 
in a bi-layer YBCO/F (a magnetic material SrRuO was 
used as the ferromagnet). They found a dip in i>{e) at en- 
ergies e = eV < IQmeV in bi-layers with the thickness of 
the ferromagnet < 26nm, whereas the "magnetic" length 

was estimated to be 3nm, that is, much shorter than 
the thickness of the F layer. Another interesting finding 
was that the change of DOS occured only in the vicin- 
ity of DW. As a possible reason for the observed effects, 
the authors of Ref. fll!| considered LRTC. However, be- 
lieving that this component had to spread all over the 
sample rather than to be located near DW, they finally 
ruled out this possibility. 

Independent of whether LRTC was observed in Ref. 
ll\ or not, it is clear that theory of the HTS/F het- 
erostructures is in a great demand. Although the super- 
conducting pairing in both conventional superconductors 
and HTS cuprates is singlet, one cannot just use for HTS 
the theory developed previously for the conventional su- 
perconductors because of a special symmetry of the order 
parameter A in these materials. 

The HTS cuprates are layered compounds with a weak 
coupling between the layers. The symmetry of the order 
parameter in HTS differs from that in the BSC supercon- 
ductors In the simplest version the order parameter 
A in HTS has the structure: A(p) = Aq{pI -pl)/p^ = 
Aocos2((/) — 4)q), where 4> is the azymuthal angle in the 
(a, b) plane normal to the c crystallographic axis, Pa and 
Pb are projections of the momentum on the axes a and b. 
Such a dependence of A(p) in HTS leads to many inter- 
esting phenomena specific to HTS. For example, the gap 
in the excitation spectrum turns to zero for certain direc- 
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tions. In addition, the sign-variable angular dependence 
of A is the reason for the sign change of the Josephson 
current in tunnel junctions composed of two HTS with 
different orientations of crystallographic axes Q . At the 
same time, we are not aware about any discussion of pe- 
cuharities of PE in HTS/F heterostructures distinguish- 
ing them from the PE in conventional S/F systems. 

In this Letter we consider proximity effects of the HTS 
in contact with disordered ferromagnets containing do- 
main walls. We concentrate on the most interesting sit- 
uation when the c-axis is perpendicular to the interface, 
such that the layers of the HTS are parallel to the inter- 
face. This geometry corresponds to the experiment 
and is special because, in cases when c— axis is not or- 
thogonal to the interface, PE is similar to the one in an 
s- wave S/N (or S/F) bi-layer [H, [3, E | . 

We will see below that in the geometry considered the 
superconducting d— wave condensate penetrates a normal 
metal with the mean free path / over a distance of order Z, 
which is much shorter than the length ^^v, characterizing 
the penetration of the s— wave superconductivity. One 
can say that PE is almost absent in this case but, as we 
will show, it is restored if one replaces the normal metal 
by a ferromagnet with domain walls perpendicular to the 
interface. In this case an odd triplet s—wave component 
is induced near DW and this is just the LRTC predicted 
in Ref. ji] for S/F heterostructures with conventional su- 
perconductors. It penetrates over much longer distances 
of order ^jv • Previously, inhomogeneities of the exchange 
field lead to generation of the odd triplet s—wave con- 
densate from the conventional s—wave singlet supercon- 
ductor. Now the presence of DW leads to formation of 
the same s—wave component but starting from the d— 
wave singlet superconductivity. 

Before we give details of our calculation, let us explain 
qualitatively why this interesting effect occurs. First, 
part of the effect, namely the destruction of the proxim- 
ity with the d —wave superconductor and normal metal is 
quite simple. It is well known the impurities destroy the 
d-wave pairing. In the case, when the c-axis is perpen- 
dicular to the interface, the symmetry of the condensate 
penetrating into the normal metal is the same as that in 
the superconductor and it is destroyed in completely the 
same way. 

Adding a domain wall perpendicular to the interface 
brings a coordinate dependence of the magnetic momen- 
tum and therefore leads to a distortion of the momentum 
dependence of the superconducting condensate. One can 
say that the symmetry of the condensate function in F 
is broken because all quantities depend on the orientaion 
of the rosette in S with respect to the plane of the DW. 
The new distorted rosette can be expanded in angular 
harmonics but now there is no symmetry reason for van- 
ishing the s-wave component. Then, as previously for 
the case of the s-superconductors with a DW wall, both 
the singlet and odd triplet component of this s-wave are 
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FIG. 1: Schematic picture of a HTS/F bi-layer with a domain 
wall DW of the Neel type (the shadowed stripe). The signs 
in the circles denote the direction of the magnetization in 
domains. 



generated and the latter can penetrate the normal metal 
over long distances along the DW. 

Our theory is based on using the method of quasiclas- 
sical Green functions. Within this approach one has to 
solve a proper equation [lH for the Green functions and 
then calculate physical quantities. 

In order to simplify the problem we assume a weak 
PE, which corresponds to a small transmission through 
the HTS/F interface. In this case the amplitude of the 
condensate function / in the ferromagnet is small and 
the Eilenberger equation can be linearized. As a result, 
this equation takes in the ferromagnet the following form 
(see for example 0,01) 



s^nw • Z(n • V)f3(g)/ + n^f - iAft,{cosa(y)[<T3, /]+ + 
+ sin«(y)f3<»[<T2,/]} = (/(n,r)), (1) 

where = 1 -I- 2|w|r, = hrsgnuj^uj = TTT(2n + 1) 
and n — p/p is the unit vector parallel to momentum. 
The products hr, \uj\t are assumed to be small. The an- 
gle brackets mean the angle average, [(73,/]+ and ['J2,/] 
denote anticommutator and commutator. 

We consider the case when the magnetization is ori- 
ented along the z— axis far away from the DW and rotates 
in the DW in the {y, z} plane changing asymptotically 
its sign (see Fig. 1). This rotation is described by the 
angle a{y). The c-axis is assumed to be orthogonal to 
the interface. 

The quasiclassical Gor'kov function / entering Eq. ([T]) 
is a 4 X 4 matrix in spin (ct^) and particle-hole (fi) 
spaces and it depends on coordinates r ={x, y} (here 
the subindex i = 0,x,y, z and ctq, tq are the unit matri- 
ces) . This matrix is off-diagonal in the particle-hole space 
and contains both the diagonal (/±) and off-diagonal el- 
ements in spin space. We use the boundary conditions 
derived in Refs. Il6l 



sgn{uj ■ n.j:){f{nx,ny) - f{-n^,ny)) = T{nx)h®fs (2) 
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where T{nx) is the transmission coefBcient, and fs = 
CF^®fs, fs = 7^2/5 is the quasiclassical Gor'kov func- 
tion in the HTS. This function is assumed to be un- 
perturbed by the weak PE and has the form: fs = 
A(n)/v/^^2 ^ A2(n) with A(n) = Aoipl - pD/p"^. Solv- 
ing Eq. II]) for / with the boundary conditions, Eqs. 
(HI), one can write the condensate contribution (5i^(e) to 
the DOS at the outer surface of the F film as Si'{e) ~ 
-(l/8)Tr(fo(^<7o(/'(n,x)))U=_,e.x=d 

The solution of Eq. (|T| can be found in a way similar 
to the one used in Ref. [6|. First, we make a transforma- 
tion introducing a new matrix / = UfjjU^, where the 
transformation matrix U = exp(if3(Tia(t/)/2) describes a 
rotation in the spin space. Then, we represent the matrix 
fu as a sum of symmetric and antisymmetric in momen- 
tum space parts: fu = § + d. Substituting the matrices 
§ and d into Eq. ([T]), one can express the antisymmetric 
matrix in terms of the symmetric one. For example, for 
the elements a± diagonal in the spin space we obtain 



a± • K/i± = —sgnuj ■ l{{n ■ V)f3(g)s± + iQuySi} (3) 

where Q = da{y)/dy, Kh± = k^j ± 2i\h- It is clear from 
Eq. ([3]) that the parity with respect to the Matsubara 
frequency lo is different for the antisymmetric part d± 
that for the symmetric ones s±,si. 

The off-diagonal component di describing LRTC can 
be found from a similar equation obtained from Eq. ([T]) 
by making the replacement Kh± — > n^j and s± —f si, 
si — > So on the right-hand side, where sq = [s+ -I- S-]/2. 

As concerns the diagonal {s±) and off-diagonal (si) 
symmetric parts, they can be found from the equations 



+iny [/2(n- V)(Qf3(g)Si 



Kh±{s±) + 
lQKh±sgnujdi], (4) 



+iny [f{n ■ V){QT3i^sn) — HQK^sgnivdo], (5) 

where oq = [0+ + a_]/2. 

Using Eqs. ^ and ^ one can write the boundary 
conditions as 



s± = ± — — jsexp — 

l\n^ 



(7) 



The matrix s± determines the amplitudes of the singlet 
and triplet, Sz = 0, components. The angular average 
(s±) is zero because (fs) — 0. The amplitude of the sin- 



glet component is proportional to S3 = (s_| 



_)/2 and 



the amplitude of the triplet component with zero projec- 
tion of the total spin on the 2;— axis is proportional to Sq. 
Both the components decay over a distance of order /. 

Obviously, this statement remains valid in the case of 
a HTS/N structure. Putting ft. = in Eqs. we see 
that only the components S3 and d^ are not equal to zero 
and that the condensate penetrating the normal metal 
preserves the d— symmetry and decays on the mean free 
path which is much shorter than the length ^^r. Of 
course, this is valid only when the axis c is perpendicular 
to the interface. Any deviation from this direction results 
in generation of the s— wave component penetrating the 
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normal metal over the distance 

Now we come to the demonstration that the presence of 
the ferromagnetism and DW restores the long range pen- 
etration of the superconducting condensate. This phe- 
nomenon occurs because such an inhomogeneous con- 
figuration of the exchange field induces the odd triplet 
s— wave component of the superconducting condensate 
with the projections of the spin Sz — ±1- 

Substituting the matrices Sq and Oq into Eq. ([5]), one 
can obtain the equation for the triplet, = ±1, compo- 
nent, si. We represent si in the form: si = siav + si^ 
with Slav = (si), such that, siav does not depend on an- 
gles and (si^) = 0. Equations describing siav and Si^ 
have quite different forms. Since (si^) = 0, one can see 
from the left-hand side of the equation for si^ that the 
characteristic scale for decay of si^^ is the mean free path 
I (in the considered dirty limit « 1). In the equation 
for Slav the second term on the left hand side and the 
first term on the right hand side almost compensate each 
other and one has: n'^siav — i^uj{si) ~ (2|a;|T)sia„- There- 
fore the angular averaged part of si decays over a large 
distance of the order of 1/ y/\uj\T ~ y^D/T and this is 
just the Sz = ±1 odd frequency triplet component. 

In the main approximation the equation for siav reads 



%x|^^U=o = T^^r(n^)/5, asi/(9a;|^=o,d = 0. (6) 
ox 2 

Eqs. completely describe the penetration of the 

superconducting condensate for the geometry in Fig.l 
that may correspond to the experiment 



We solve 



these equations assuming that IQ <^ 1, which means that 
the mean free path I is much smaller than the width of 
the DW, w (Q w~^). Then, in the main approximation 
we find from Eqs. ^ and ([5]). 



I 



■ d'^siaviq) 



+ K„siav{q) = 'iil {dQ/dy)q{n sq{x)), 



(8) 

siav{q) = J dysiaviy) exp{iqy) and = 6\uj\T + q'^P. 

The function sq {x) decays over a distance of the order 
I, whereas the characteristic scale for siavix,y) is much 
longer: 



l/\/6ujT « y/D/QjrT. Therefore, the term on 
the right hand side can be replaced by BS (x) , where B 
does not depend on x. Solving Eq. ([5]) in this approxima- 
tion and Fourier transforming back to y, one can obtain 
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correction to the DOS 




FIG. 2: a) Normalized correction to the DOS vs normalized 
energy e/Ao for d = 2^A,m = AoTm = 0.3 (solid curve) and 
d = ^AyTn = 0.6 (point curve), where — _D/(2Ao); b) De- 
pendence of LRTC amplitude on the normalized y-coordinate 
at X = d for w = ^ui,d = 0.5£,ui (point curve); w = ^:j,d — 
(upper solid curve) and w = 0.2^u, d — (lower solid curve). 
The y-coordinate is measured in units ^ui- 



the function siav{d,y) that determines the contribution 
5v{e) of the superconducting condensate to the density of 
states on the outer side of the ferromagnet. The solution 
for the LRTC siav{x^ y) at the x = d can be presented in 
the form 



Note that spin-dependent scattering makes the char- 
acteristic length of the LRTC decay shorter @] . The 
amplitude of LRTC contains a small parameter {p /wd). 
If this parameter is not small, the formulas given above 
are valid qualitatively. 

We believe that our results can help to explain the 
experimental data on the density of states in a setup 
analogous to that considered here ll|, but more infor- 



mation about the DW, barrier transparency T(nx) etc., 
is needed for a detailed comparison. 

To conclude, we have demonstrated that the presence 
of a domain wall at the interface of a high tempera- 
ture, d— wave superconductor/ferromagnet bi-layer leads 
to generation of an odd frequency triplet s— wave compo- 
nent of the condensate. As a result, the superconducting 
condensate penetrates the ferromagnet along DW over 
distances much exceeding the penetration depth of the 
(i— wave condensate into a normal metal. This explicitly 
demonstrates enhancement of the proximity effect by an 
inhomogeneous exchange field. 

We thank SFB 491 for financial support. 



Slav{d,y) = {i/ ^)ifi\hC0s{2ifQ)TavfSavS{y) (9) 

where Tav and fsav arc the transmission coefficient and 
condensate function in S averaged over angles with corre- 
sponding weights; the Fourier component of the function 
S{y) is equal to: S{q) = P[dQ/dy)q[KqS\n\i{nqd/l)\^^ . 

We emphasize that, in contrast to the singlet compo- 
nent /5, the LRTC, siau, is an odd function of the Mat- 
subara frequency because \h ^ sgnuj (see Eq.(IT])). Thus, 
the odd-frequency triplet s— wave component siav arises 
as a result of the influence of the non-homogenous mag- 
netization in the DW on the c?— wave singlet component. 
One can see that the LRTC is zero if the nodes in the 
spectrum of the d— wave superconductor lie in the plane 
of the DW. 

One can express the spatial dependence of this com- 
ponent in an explicit form in limiting cases. If DW is 
broader than the LRTC penetration depth we can write 
for the correction 6i'{e) — i'(e) — 1 to the density of 
states: Sueffie) = Siyie)[i3/2)PXhTavidQ/dy)]-\ We 
calculated the function v^f f (e) numerically and repre- 
sented in Fig. 2a an effective local correction to the DOS 
6vf,ff{e) aX X — d caused by the LRTC. 

In the opposite limiting case a simple form for the 
LRTC can be obtained if the DW is approximated by 
a step-like function {Q = n/w for \y\ < w/2 and Q = 
for \y\ > w/2). In this approximation we calculate the 
spatial dependence of Seffiy) — —{w/l)S{y), for various 
values of w and d (Fig. 2b). We see that the LRTC decays 
exponentially from DW over a long distance of order ^jv. 
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